Shortcuts to adiabaticity by superadiabatic iterations 
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Different techniques to speed up quantum adiabatic processes are currently being explored for 
applications in atomic, molecular and optical physics, such as transport, cooling and expansions, 
wavepacket splitting, or internal state control. Here we examine the capabilities of superadiabatic 
iterations to produce a sequence of shortcuts to adiabaticity. The general formalism is worked out 
as well as examples for population inversion in a two-level system. 
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I. INTRODUCTION 

There is currently much interest to speed up quantum 
adiabatic processes for applications such as fast cold- 
atom or ion transport, expansions, wave-packet split- 
ting or internal state population and state control 
Different techniques have been put forward and/or ap- 
plied. Among them Dcmirplack and Rice [2|-y|, and 
Berry [j| proposed the addition of a suitable counter- 

diabatic term H^)(t) to the time dependent Hamilto- 
nian Ho(t) whose adiabatic dynamics is to be imple- 
mented. With that term added transitions in the in- 
stantaneous cigenbasis {\no(t))} of Ho(t) are suppressed, 

H (t)\n (t)} = E^\t)\n (t)), while there are in general 
transitions in the instantaneous eigenbasis of the new full 
Hamiltonian Hq + . Experiments that implement 
these ideas have been recently performed in different two- 
level systems @, 0] • The same HfJ (t) also appears nat- 
urally when studying the adiabatic approximation of the 
reference system, the one that evolves with Ho(t), see e.g. 
[H . The reference system behaves adiabatically, following 
the eigenstates of Ho(t), when the countcrdiabatic term 
is negligible, and the adiabatic approximation is close to 
the actual dynamics. This is made evident in an interac- 
tion picture (IP) based on the unitary transformation 
A (t) = J2 n \no(t))(n (0)\. (The "parallel-transport" 
condition (no(t)\ho(t)) = is assumed hereafter to define 
the phases.) From the Schrodingcr equation ihdtipo(t) 
II ' and defining -01 (t) = A^ipo^), the IP equa- 

tion ihd t tpi(t) = Hi(t)tl)\{t) is deduced, where H\{t) = 
A\{t)(H (t) -K (t))A {t) is the effective IP Hamiltonian 
and K (t) = ihA (t)Al(t) is a coupling term. If Ko(t) 
is zero or negligible, Hi(t) becomes diagonal in the basis 
{|rio(0))}, so that the IP equation becomes an uncoupled 
system with solutions 



hM*)> = iMt)h&i(o)>, 



where 



0i(t) = £)|no(O))e 



(1) 



S(0)| (2) 



is the unitary evolution operator for the uncoupled sys- 
tem. Correspondingly, from |VoC*)) = -^oWl^iM), 

=El»o(t))e-^»' B « )(t ' )dt '(n„(0)|*(0)), (3) 



where we have used \ipi(0)} = \ipo(0)) since A\(0) = 1 
by construction. The same solution, which, for a non- 
zero Ko(t), is only approximate, may become exact by 
adding to the IP Hamiltonian the countcrdiabatic term 
AQ(t)Ko(t)Ao(t). This requires an external intervention 
and changes the physics of the original system, so that 
Ui (i) describes the evolution exactly. In the IP the modi- 
fied Hamiltonian is = H^t) + Al(t)K a (t)A (t) = 
Al(t)H (t)A (t), and in the Schrodingcr picture (SP) 
the additional term becomes simply Kq. The modified 
Schrodinger Hamiltonian is Hg(t) = Ho(t) + Ko(t), so 

we identify H^(t) = Ko(t). A "small" coupling term 
Ko(t) that makes the adiabatic approximation a good 
one also implies a small counterdiabatic manipulation 
but, irrespective of the size of K (t), H { 1] (t) provides 
a shortcut to slow adiabatic following because it keeps 
the populations in the instantaneous basis of Ho(t) in- 
variant, in particular at the final time tf. Moreover, if 

K (0) = and K Q {t f ) = 0, then H { 1] = H at t = 
and t = tf. This is useful in practice to ensure the con- 
tinuity of the Hamiltonian at the boundary times: usu- 
ally H (t < 0) = F (0) and H (t > t f ) = H (t f ), so 
Ko(t < 0) = Ko(t > tf) = 0, i.e., H Q (t) is the actual 
Hamiltonian before and after the process. 

The previous formal framework may be repeated itera- 
tively to define further IPs by diagonalizing the effective 
Hamiltonians of each IP. These iterations were used to 
establish generalized adiabatic invariants and adiabatic 
invariants of n-th order by Garrido 0. Berry also used 
this iterative procedure to calculate a sequence of cor- 
rections to Berry's phase for cyclic processes with finite 
slowness, and introduced the concept of "superadiabatic- 
ity" 11011- For later developments and applications see e.g. 

a EMU- 

The idea of superadiabatic iterations is best under- 
stood by working out explicitly the next interaction pic- 



2 



turc: 1 let us start with iHdtipi{t) = Hi(t)ipi(t) and treat 
it as if it were, formally, a Schrodinger equation. The di- 
agonalization of Hi(t) provides the eigenbasis {\ni(t))}, 

Hi(t)\ni(t)) = En\t)\n\(t)) , that we fix again with the 
parallel transport condition, (ni(t)\hi(t)) = 0. A new 
unitary operator A\ = J2 n l n iW)( n i(0)l plays now the 
same role as Aq in the first (adiabatic) IP. It defines a 
new interaction picture wave function ip 2 (t) = ^iW^iW 
that satisfies ihd t ^p2{t) = H 2 (t)tp 2 (t), where H 2 (t) = 
A\(t)(H-L(t) - Ai(t))Ai(i) and R\ = ihAiA[. If K x is 
zero or "small" enough, i.e., if a (first order) superadi- 
abatic approximation is valid, the dynamics would be 
uncoupled in the new interaction picture, namely, 

\i> 2 (t)) = u 2 (t)\MQ)), (4) 

where 

U 2 (t) =^|n 1 (0)) e -^^ £ " , ( t '^'(n 1 (0)| (5) 

n 

is the approximate evolution operator in the second IP 
for uncoupled motion. It may happen that a process 
is not adiabatic, since K$(t) may not be neglected, but 
(first-order) superadiabatic when K\ (t) can be neglected. 

(2) 

Transforming back to the Schrodinger picture, (t)) = 
A)(t)Ai(i)t/ 2 (i)IV'2(0)) becomes 

l4 2) W) = EEl m oW)( m o(0)|m(t)) e -^o £ " )(t ' )d *' 

n m 

x (ni(O)l^o(O)), (6) 

and |^ 2 (0)) = |Vo(0)) since A o (0) = A^O) = 1. Garrido 
distinguished two different aspects 

• Generalized adiabaticity: The evolution operator 
Ao(t)Ai(t)U 2 (t) provides an approximation to the 
actual (Schrodinger) dynamics up to a correction 
term of order 1/tj. This is so without imposing 
any boundary conditions (BCs) at t = and t = tt 
on the Hamiltonian Hq. 

• Higher order adiabaticity: Ao(t) A\(t)U2(t) does 
not guarantee in general that |rio(0)) evolves into 
\no(tf)), up to a phase factor. If this is the objec- 
tive, in other words, if a superadiabatic approxima- 
tion should behave, at final times, like the adiabatic 
approximation, up to phase factors, then some BCs 
have to be imposed. Garrido discussed how gener- 
alized adiabaticity implies higher order adiabaticity 
when BCs at the boundary times are imposed on 
the derivatives of Hq. 



1 The first IP and iteration just described, with dynamics governed 
by Hi(t), generates the modified dynamics based on Hq in the 
SP. This iteration may be naturally termed as "adiabatic" , since 
the unitary transformation used, Ao, relies on the usual adiabatic 
basis. Moreover this is the IP used to perform the adiabatic 
approximation by neglecting Kq. The second iteration may be 
considered as the first "superadiabatic" one. 



(Garrido's distinction does not apply in since there, 
it is assummed from the start that all the derivatives of 
the Hamiltonian Hq vanish at the (infinite) time edges.) 
The second aspect is crucial to design shortcuts to adia- 
baticity for finite process times using the superadiabatic 
iterative structure, so let us be more specific. First no- 
tice that Eq. ([5]) becomes exact if the term A\K\A\ is 
added to the IP Hamiltonian, so that now the modified 
IP Hamiltonian is = H 2 + A\KiAi = AIh^. 

Then the modified SP Hamiltonian becomes Hq (t) = 

Ho(t) + H^ d \t), where H$(t) = A (t)Ki(t)Al(t). How- 
ever, quite generally the populations of the final state © 
in the adiabatic basis {\n (tf))} will be different from 
the ones of the adiabatic process, unless (a) {|no(0))} = 
{\ni(tf))} 7 up to phase factors, and (b) {|«o(0)}} = 
{|ni(0)}}, also up to phase factors, (a) is satisfied when 
Ko(tf) = 0. This makes H\(tf) diagonal in the basis 
{K(0)}} and [#!(£/), iJ o (0)] = 0. (b) is satisfied when 
Ao(0) = 0, which implies iii(O) = iAi(O). In summary 
the requirement is that the eigenstates of Hi(t) at t = 
and t — tf coincide with the eigenstates of Hq(0). If, 
in addition, A'i(0) = K\(tf) = 0, not only the popula- 
tions, but also the initial and final Hamiltonians are the 
same for the "corrected" and for the reference processes, 
namely ^ 2) (0) = H (0) and H^\t f ) = H {t f ). 

Further iterations define higher order superadiabatic 
frames with IP equations ihd t ^j{t) = Hjijjj(t), where 

Hj = ^j-i(Hj-i - Kj-i)Aj-i, (7) 
Kj = ihAjA], (8) 

with A 3 {t) = £,>*(*)) <»j (0)1 and H 3 {t)\n 3 {t)) = 
En\t)\n 3 {t)) . As the A 3 (0) = 1 by construction, there 
is a common initial state \i[>j(0)) = |V'o(0)) for all itera- 
tions. The general form for the modified IP Hamiltonians 
is ff« = Hj + .I] ,/\, i = Al^Hj^Aj^!. Thus, 

the form of the modified Hamiltonians in the SP is 

H^^H. + H^, (9) 

where the SP counterdiabatic term is 

//,.; n.hjf ihnJ.An:. (10) 

with B 3 = A ■ ■ ■ Aj_i and B = I. If K 3 -i(t) is 
small or negligible, H 3 (t) becomes diagonal in the basis 
{|nj_i(0))} and the IP equation becomes an uncoupled 
system with solutions \ip 3 {t)) = U 3 (t)\ipj(Q)}, where 

C/ 3 (t)^| Vl (0))e-i/« E ^ 1) ( t ')' ,t '( rij _ 1 (0)| (II) 

n 

is the approximate evolution operator in the j-th IP. Cor- 
respondingly, the approximate solution in the SP is given 

by \^\t)) = Ao(t)Mt) ■ ■ ■ Aj-i(t)Uj(t)\M0))- This 
solution becomes exact if the H^ 1 ^ term is added to 
Hq, where in general the populations of (tf)) in the 
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adiabatic basis {|fto (*/))} wn l be different from the ones 
of the adiabatic process, unless appropriate BCs are im- 
possed. These boundary conditions are made explicit in 
the next section, and correspond partially to the condi- 
tions discussed by Garrido in Q to define "higher order 
adiabaticity" . 

Is there any advantage in using one or another counter- 
diabatic scheme? There are several reasons that could 
make higher order schemes attractive in practice: one 
is that the structure of the H^) may change with 
j. For example, for a two- level atom population in- 
version problem, H^J = h(Q /2)a y , whereas H^J = 

7i(0i/2)(cos <do<r x — sin 0oo\z), where the Qj are the polar 
angles corresponding to the Cartesian components of the 
Hamiltonian Hj, and the a Ul with u = x,y,z, are Pauli 
matrices 0. (We shall use the Cartesian decomposition 
Xa x + Ytjy + Za z for different Hamiltonians below.) 

A second reason is that, for a fixed process time, the 
cd-terms tend to be smaller in norm as j increases, up 
to a value in which they begin to grow (l5| . An opti- 
mal iteration may thus be set [TJ El ■ The "asymptotic 
character" of the superadiabatic coupling terms and the 
eventual divergence of the sequence can be traced back 
to the existence of non-adiabatic transitions, even if they 
are small [Ioj |. 

To generate shortcuts one should pay attention though 
not only to the size of the cd-terms but also to the feasi- 
bility or approximate fulfillment of the required BCs at 
the boundary times. Thus, it may happen that an "op- 
timal iteration" , of minimal norm for the cd-term, fails 
to provide a shortcut because of the BCs, as illustrated 
below in Sect. IV. 



II. BOUNDARY CONDITIONS FOR 
SHORTCUTS TO ADIABATICITY VIA 
SUPERADIABATIC ITERATIONS 

In this section we set the boundary conditions that 
guarantee that H^\t) provides a shortcut to adiabatic- 
ity. We have seen that for j = 1 no conditions are re- 
quired. For j — 2 we need that = {|«o(0))} 
and {|ni(0))} = {|no(0)}}, (as before in these and similar 
expressions in brackets, the equalities should be under- 
stood up to phase factors), i.e., K {tf) = K {Q) = 0. 
For the iterations j > 2 we need that (a) {|nj_i(0)}} = 
{|n (0))}, which occurs when Kj-2(0) = Kjsify = ••• = 
Kx(0) = K (0) = 0, and (b) {^(tf))} = {[^- 2 (0)», 

{K_ 2 (i/)>} = {k- 3 (o)», {K- 3 (t/)>} - {k- 4 (o)>}, 

... , and {\ni(tf))} = {|no(0))}. This amounts to impos- 
ing Kj- 2 (t f ) = K^ 3 (t f ) = ... = R\(t f ) = K (t f ) = 0. 
The vanishing of Kj>(0) for j' < j — 2 implies that 
ff o (0) = -Hi(O) = ... = #j_i(0), so (a) and (b) com- 
bined may be summarized as {|?v(0))} = {\ n j'(tf))} = 
{|no(0))} for all f < j — 1. Garrido showed that cancel- 
ing out the first Z-th time derivatives of Ho(0) and Ho(tf) 
makes Kj(0) = and Kj(tf) = 0, for j = 1, / — 1, re- 
spectively [9(. However canceling out the derivatives of 



Hq is a sufficient but not a necessary condition to cancell 
the coupling terms, so we find it more useful to focus in- 
stead on the coincidence of the bases, this is exemplified 
in Sect. IV. 



III. ALTERNATIVE FRAMEWORK WITH A 
CONSTANT BASIS 

An alternative to the formal framework described so 
far provides computational advantages. It was implicity 
applied by Demirplak and Rice for a two- level system [4| . 
We shall here generalize and formulate explicitly this ap- 
proach and show its essential equivalence to the former. 
The main idea is to use instead of the Aj a different set 
of unitary operators, Aj(t) — ^2 n \hj(t))(n\, to define 
the sequence of interaction pictures, where \hj(t)) are 
cigenstates of the new IP Hamiltonians Hj (t) , such that 
Hj(t)\hj(t)) = En (t)\fij{t)), and {|n)} is a constant or- 
thonormal basis equal for all j, which in principle does 
not necessarily coincide with \rij(0)). Similarly to Eq. 

H^A^H^-k^A^t, (12) 

where Kj = iHAjAj = ihJ2n The coun- 

terdiabatic terms in the SP are introduced as before, 
= BjKjBt, where Bj = A ■ ■ ■ A,_i with B Q = I. 
We shall next show that these cd-terms are independent 
of the chosen constant basis, so that H^)(t) — H^)(t). 
Therefore, it is worth using Aj (t) instead of Aj (t) since 
they are simpler operators and significantly facilitate the 
manipulations as a common basis is used. 

Let us start with the first iteration. Since Ho(t) = 
H (t), then Ei°\t) = E [ °\t), \h (t)) = |n (t)), and 
A'o = K Q , so H { c ° d ] = Hj$. In addition, from Eq. Q, 
Hq — Kq = Ao-HiAq, and substituting it in Eq. (|T2"j) leads 
to 

Hi = uoHml (13) 
where we have defined a constant unitary operator 
m = AlA = J2\n}(n (0)\, 

n 

iio = 0. 

Using 

Hj(t) =J2Mt))El j) (t)(n j (t)\ (14) 

n 

and 

^•W = ^|n J -(t))^(t)(n j (t)|, (15) 

n 
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for j = 1 in Eq. JT3), wc get that E^ (t) = E^\t) and 
\ni(t)) = M |ni(t)), while \n) = w o |n o (0)). Expanding 
H^y = AqKxAL we have that 



cd 
AD 



H$(t) = ih 



J2 I "o (*)> (n |^i (t)> (m|/) (fi (t) |p) (po (t) | . 



Using now (m|Z) = <5 T „z, |no(t)) = \no(t)), \n) = 
Mo|no(0)}, and \ni(t)) = uo\ni(t)), it follows that H^J = 
H%. Also, Ki = AlAoK^Ao = u K lU l 

Repeating these steps for j ^ 1, Hj = Uj-iHjU^_ 1 
and Kj = Uj_iKjiij_ 1 , where 

Uj = A^Uj-iAj = ^2 \ n ) i n 3 (°) I . 

n 

Uj = 0. 

This leads to E#\t) = E [ J\t), \hj(t)) = u j _ 1 |n J -(t)) ) 
and |n) = «j_x|n./_i(0)}. Thus, for all j > 0, 

H cd - H cd ■ 

The boundary conditions to achieve shortcuts to adia- 
baticity take the same form as for the original frame- 
work in the previous section. Since Kq = Kq and 
Kj = Uj-iKjUj_i for j 1, for the j-th iteration, with 
j > 1, we need that K (0) = l<i{0) = ... = Kj- 2 (0) = 0, 
and Ko(tf) = Ki(tt) = ... = Kj^itf) = 0. Let us recall 
that no conditions were required for j = 1, although, as 
shown in the next section, using a convenient (constant 
or initial adiabatic) basis for specific Hamiltonians may 
also lead to conditions for j = 1. 



IV. TWO-LEVEL ATOM 

The general formalism will now be applied to the two- 
level atom. Assuming a semiclassical interaction between 
a laser electric field and the atom, the electric dipole and 
the rotating wave approximations, the Hamiltonian of 
the system in a laser-adapted IP (that plays the role of 
the Schrddinger picture of the previous section) is 



H (t) 



h ( -A(t) n R (t) 

fl R (t) A(t) 



(16) 



where fl R (t) is the Rabi frequency, assumed real, and 
A(£) is the detuning, in the "bare basis" of the two level 
system, |1) = (I), 2) = ("). The Hamiltonians of the 
consecutive interaction pictures can be written as j!| 



Hj(t) 



Zj(t) Xj{t)-iYj{t) 
Xj{t)+iY 3 (t) -Zj(t) 



(17) 



or Hi 



Xj(7 x 



ZjO z 



[19|. Then, X (t) = 
Hfl R (t)/2, Y (t) = and Z (t) = -HA(t)/2. Xj, Yj, 
and Zj are the Cartesian coordinates of the "trajectory" 



of Hj{t). It is also useful to consider the corresponding 
polar, azimuthal and radial spherical coordinates, @j(t), 
$j(t), and Rj(t) HE1, that satisfy 

COS ( 0J .) = |l, sin(9,) = < Qj < 7T, 



R, 



= ^, Biii(* i ) = 5-. 0<$ i <27r, (18) 



with Rj = JXj + Yf + Z] and Pj = JXj + Yj, where 
the positive branch is taken. The eigenvalues of Hj(t) 

are E{ = —Rj and E^ = Rj, and the corresponding 
eigenstates {\fij(t))} are 



|2;} 



e, 



|l)-e 



i<S> j/2 



I ^ ) 



f ^)|l)+e^/ 2 sin^ ) 



(19) 



where the phase 
e,-(t) = 



1 

2 7o 



^ (t') cos [9j ;(*')]<&' (20) 



is introduced to fulfill the parallel transport condition 
(hj\hj) = 0. Wc define Aj = \lj(t))(l\ + [ 2_y (*) > <2 1 . The 
matrix Aj (t) under these conditions is 



At 



sm |& 1^-^/2 cos f ®L) e -iei-i^/2 



COS 



^ L \ p -re j +i4> j /2 



Then, from Eq. 



(21) 



— Qj sin ^ cos sin (26j 

6j cos - sin ($j) sin (29j) 



~2 



^sin 2 (9,-)^. 



(22) 



Note that A^-K^A; = A'-Aj has only non-diagonal ele- 

f — L 

ments in the bare basis {|1), |2)} [J]. 

From Eq. ||7J), the Cartesian coordinates of Hj +1 (t) 
are 





?! 




A J+ i - 


2 


9 




?! 




= 


2 





= 





-9j cos (2ej) - 4>j sin (9j) sin (2ej) 



(23) 



In general, if $j (t) is constant for a particular j — J, then 
<j>j(t) = 0, and from Eq. ([2"0]). ej(t) = 0. Thus, taking 
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into account Eq. (f23|) . we have that Xj + i(t) = and 
Y J+1 (t) = -HQj/2. Equation (HHJ) leads to $j + i(t) = 
{7r/2,37r/2}, with tt/2 when Y J+1 > (Qj < 0), and 
3tt/2 when Y J+1 < (Qj > 0). If Y J+1 = 0, $ J+ i 
is discontinuous, and 9,/+i = it. Therefore, e,j+\(t) = 
{0, ±7r/2}. From here, several general conditions can be 
deduced for j' > J: <Iy>,/(i) = {tt/2, 3tt/2}, ej>>j(t) = 
{0,±^/2}, -X>>j(t) - 0, and y J+x (t) = -HQj/2 or 
fy>J+i(*) = ±h®j>-i/2. Moreover, from Eq. (1221) . 
Kj'>J = ±(M) 3 ' /2)a x with positive sign if $ 3 / (i) = 37r/2 
and negative sign if $j'(f) = 7r/2. Eq. ([TU) an d Y (t) = 
imply $„(*) = if X (i) > and $ (*) = tt if X (i) < 0. 
We may thus take J = and apply the above relations, 
for example <&o(i) = and £o(t) = 0. 2 

As we mentioned before, the method fails as a shortcut 
to adiabaticity when the boundary conditions are not 
well fulfilled. In order to have a shortcut generated by 
the iteration j we require that A(i) and tlji(t) are such 
that 
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(24) 
(25) 



for < j' < j, up to phase factors. For j' = a natural 
and simple assumption is that the bare basis coincides 
initially with the adiabatic basis, i.e., Eq. (f2~4")k at tf 
we assume that the bare and adiabatic bases also coin- 
cide, allowing for permutations in the indices and phase 
factors. 

At t = 0, using Eq. (fT9"j) . taking into account that, 
from Eq. (gDJ), £,'(0) = 0, and that $ (*) = 0, 
sin[6 3 ,(0)/2] = 1 and cos [8 3 > (0)/2] = arc required, 
or 9 3 v(0) = 7T. Then, cos[0 3 /(O)] = Z r (0)/ R r (0) = -1. 
This condition is fulfilled if 
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FIG. 1: (Color online) The (a) X, (b) Y and (c) Z components 
for the Landau-Zener scheme, of: Ho(t) (black dots), Hq 



(purple solid line), -ffp 2 ' (orange dotted line), Hq 31 (red dashed 



(3) 



4(o)»x 2 (o) + y 2 (o), 



(26) 



line), and Hq (blue dot-dashed line). In (a) and (c) the black 
dots and the green thin line coincide and in (b) the black dots 
coincide with the orange dotted line. Parameters: a = —20 



as long as Z> =0 (0) < 0, and knowing that Zf >0 (t) = MHz , n ,iz = 0.2 MHz, and t f = 0.2 /is. 
— Rj/-i(t) < 0. The condition (|2"6")l can be simplified for 
specific j'-values as 



At t = t 



|Z (0)| » |X (0)|, 
l^-'>o(0)| » |^>o(0)|. 



Zf,(t f )>XUtf) + YP(t f ) 



(27) 
(28) 



(29) 



should be satisfied, where now, 0rj(i/) can be either 0, if 
Z (tf) > 0, or 7r if Z (tf) < 0, and 3 '>o(*f) = tt. As 
before, this condition splits into 



As an example we consider now a Landau-Zener 
scheme for Hq (for the Allen-Eberly scheme we have 

found similar results), and study the behavior of 
with j = 1, 2, 3, 4, and the populations of the bare states 
driven by these Hamiltonians. 

For the Landau-Zener model A(t) is linear in time and 
f2fl(t) is constant, 



Aj,(t) 
fin,i*(t) 



a(t-t f /2), 
"OJzi 



(32) 



\Z (t f )\ » \X Q (tf)\, 
\Zj>>o(tf)\ » \Y r>0 (t f )\ 



(30) 
(31) 



where a is the chirp, and f2o,(z is a constant Rabi fre- 
quency. Condition ([27|) can be restated as 



t / »2 



2 The analysis in this paragraph follows closely Q], but some of 
the results differ, in particular the values allowed for the phases 
£j'>j- 



fi . 



Iz 



(33) 



We consider the parameters a = —20 MHz 2 , f2o,iz = 0.2 
MHz, and tf = 0.2 fis for which the dynamics with Hq 
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TABLE I: Maxima of the X and Y components of Ho and Hq 
for j = 1,2,3,4,5. Parameters: a = -20 MHz 2 , Q ,i* = 0.2 
MHz, and tf =0.2 fxs. 



Hamiltonian \X max \/h (MHz) \Y max \/h (MHz) 



H /h 


0.1 





< ] /h 


0.1 


49.9 




10 





H*>/h 


S.4 


2.8 


H^/K 


46.8 


28.1 




56.2 


62.8 



is non-adiabatic, see the Appendix A. Fig. [T] shows X, 

(i) 

Y and Z components of Hq and Hq, with j — 1, 2, 3, 4. 

In Figs. QJ, andQJi and in table U we see that Hq (cor- 
responding to the first superadiabatic iteration) is op- 
timal with respect to applied intensities. Moreover it 
cancells the K-component completely, which is a sim- 
plifying practical advantage in some realizations of the 
two-level system 0, [h|. From the second superadia- 
batic iteration both intensities start to increase again. 
For the parameters above, condition (|33[) is satisfied 
since tf = 20 x \£loj z /a\, but not so condition (|28j) . 

Fig. [1] shows the disagreement between Hq^ and Hq, 
at t = and t — tf for j > 1. Fig. [2] shows that only 

= H Q + H^J inverts the population of |1), Pi(t), 
whereas the rest of the Hamiltonians fail to do so. 




0.1 0.2 

t(fXS) 



FIG. 2: (Color online) Population of the state |1), Pi(t), for 
the Hamiltonians Ho(t) (black solid line with dots), Hg 1 ^ (pur- 

(2) (1) 

pie solid line), Hq (orange dotted line), Hq (red dashed 
line), and H^ (blue dot-dashed line), with the Landau-Zener 
scheme. Parameters as in Fig. CO a = -20 MHz 2 , n 0>lz = 0.2 
MHz, and tf = 0.2 fis. 



V. DISCUSSION 



In this paper we have investigated the use of quantum 
superadiabatic iterations (a non-convergent sequence of 
nested interaction pictures) to produce shortcuts to adi- 
abaticity. Each superadiabatic iteration may be used in 
two ways: (i) to generate a superadiabatic approxima- 
tion to the dynamics, or (ii) to generate a counterdia- 
batic term that, when added to the original Hamiltonian, 
makes the approximate dynamics exact. The second ap- 
proach, however, does not automatically generate short- 
cuts to adiabaticity, namely, a Hamiltonian that produces 
in a finite time the same final populations than the adia- 
batic dynamics. The boundary conditions needed for the 
second approach to generate a shortcut have been spelled 
out. This work is parallel to the investigation by Garrido 
to establish conditions so that the approach (i) provides 
an adiabatic-like approximation Q. We have also de- 
scribed an alternative framework to the usual set of su- 
peradiabatic equations which offers some computational 
advantages, and have applied the general formalism to 
the particular case of a two-level system. An optimal 
superadiabatic iteration with respect to the norm of the 
countcrdiabatic term, is not necessarily the best short- 
cut, or in fact a shortcut at all, because of the possible 
failure of the boundary conditions. 

We end by mentioning further questions worth inves- 
tigating on the superadiabatic framework as a shortcut- 
to-adiabaticity generator. For example, other operations 
different from the population control of two-level systems 
(such as transport or expansions of cold atoms) have to 
be studied. Unitary transformations may be also ap- 
plied to simplify the Hamiltonian structure making use 
of symmetries [19j . They have been discussed before as 
a way to modify the first (adiabatic) iteration [TTI. Il9j|. 
and applied to perform a fast population inversion of a 
condensate in the bands of an optical lattice @, but a 
systematic application and study, e.g. of the order with 
respect to the small (slowness) parameter, in particu- 
lar for higher superadiabatic iterations, are still pending. 
A comparison with other methods to get shortcuts, at 
formal and practical levels would be useful too. A pre- 
liminary step in this direction, relating and comparing 
the invariant-based inverse engineering approach to the 
counterdiabatic approach of the first (adiabatic) iteration 
was presented in |20j, see also the Appendix B. Finally, 
comparisons among superadiabatic iterations themselves 
have to be performed, in particular regarding practical 
aspects such as the transient excitations involved (2lj . 

We are grateful to O. Morsch and M. Berry for discus- 
sions. We acknowledge funding by Projects No. IT472- 
10, No. FIS2009-12773-C02-01, and the UPV/EHU pro- 
gram UFI 11/55. S. I. acknowledges Basque Government 
Grant No. BFI09.39. X. C. thanks the National Natural 
Science Foundation of China (Grant No. 61176118) and 
Grant No. 12QH1400800. 
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Appendix A: Adiabaticity and boundary conditions 
for the Landau-Zener protocol 

The adiabaticity condition for a two-level atom driven 
by the Hamiltonian ([To]) is [22j 



i|n a (t)| « |n(t)| 



(Al) 



where = [n R (t)A(t) - tl R (t)A(t)]/n 2 (t) and ft = 

1/ A 2 + fl R . For the Landau-Zener scheme this condition 
takes the form 



(A2) 



The inequalities that a must satisfy so that the system 
is adiabatic and also fulfills the boundary condition (|33j) 
are 



2|fi ,u|/t/ < |a| <2ft 2 



lz- 



(A3) 



Fig. [3] shows the (shaded area) region for which a 
satisfies 20|f2 ,iz|/*/ < H < 0.2ft 2 , lz when tj = 2 
/is. No such area exists in the depicted domain for 
tf = 0.2 fis. For this shorter time the critical point where 
l/tf = fto./2/lOO corresponds to fto,z 2 = 500 MHz and 
detunings of up to 5 GHz. Both may be problematic, 
as very large laser intensities and detunings could excite 
other transitions. 



4000 

X 




100 
Q s (MHz) 



200 



FIG. 3: (Color online) 0.2flo j,, (red dashed line) and 
20\Qo, iz\/tf (blue dashed line for tf = 2 [is and black dot- 
dashed line for tf — 0.2 /is). The shaded (green) area corre- 
sponds to values of a satisfying 20jf2o,iz|A/ < M < 0.2f2§ !zj 
namely, the process is adiabatic and the eigenstates at the 
boundary times are essentially the bare states. No such area 
exists for tf = 0.2 ps in the domain shown. 



Appendix B: Invariants 

The supcradiabatic sequence may be pictured as an 
attempt to find a higher order frame for which a cou- 
pling term Kj is zero in the dynamical equation so that 
there are no transitions in some basis. This would mean 
that the states that the system follows exactly have been 
found, in other words, the eigenvectors of a dynamical 



invariant I(t) [2fJ, |23|, |24| . When counter-diabatic terms 
arc added, it is easy to construct invariants for 
from the instantaneous eigenstates of -ffo(i)- However, 
quite generally this is not enough to generate a shortcut 
to adiabaticity because the boundary conditions to per- 
form a quasi-adiabatic process (one that ends up with 
the same populations than the adiabatic one) may not 
be satisfied. A way out is to design the invariant first, 
and then H(t) from it, satisfying the boundary condi- 
tions [I(t),H(t)] = at t = and t = t f, and such that 
H(0) = #o(0) and H(t f ) = H (t f ) [2,[l|. For the gen- 




FIG. 4: (Color online) X (red solid line) and Z (blue dashed 
line) components of the Hamiltonian obtained using the 
invariant-based inverse engineering method, tf = 0.2 fj,s. 



N 
« 



150 




-75 



■150. 



FIG. 5: (Color online) The components Xo(t) (red solid line) 
and Zo(t) (blue dashed line) of Ho(t), the Y(t) component of 
Hq 1 ^ (t) (orange dot-dashed line), and the X (t) (purple dotted 
line) and Z(t) (black solid line with dots) components of Hg 2 \ 
for the Landau-Zener scheme. Parameters: a — —2800 MHz 2 , 
fl 0Jz = 30 MHz, and t f = 0.2 (JS. 



cral Hamiltonian in Eq. (|I6p . a dynamical invariant of 
the corresponding Schrodingcr equation may be param- 
eterized as I20| 



I(t) 



cos7(i) sin7(t)e i/3(t) 
sin7(£)e- l/3 W -cos 7(4) 



(Bl) 



where v is an arbitrary constant with units of frequency 
to keep I(t) with dimensions of energy From the invari- 
ance condition for /, 



dl(i) _ dl(t) 



dt 



dt 



:[I(t),H (t)} = 0, 



(B2) 
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0.1 0.2 

tfjis) 

FIG. 6: (Color online) Population of |1>, Pi(t), for the Hamil- 
tonians Ho(t) (red dashed line), Hq (blue solid line), and 
Hq 2 ^ (black dot-dashed line), with the Landau-Zener scheme. 
Parameters as in Fig. [5] a = -2800 MHz 2 , Q ,iz = 30 MHz, 
and tf = 0.2 /xs. 



the functions j(t) and j3(i) must satisfy the differential 
equations 

7 = n R sm/3, 

P = A + f7 H cos/3cot7. (B3) 

To achieve a population inversion, the boundary values 
for 7 should be 7(0) = and j(tf) = ir. Assuming 
a polynomial ansatz [20|, HH, [26| for -f(t) and f3(t), as 
7(t) = X)n=o a ni™ with the boundary conditions 7(0) = 



tt, 7 (^) = 7(0) = 7(*/) = 0, and p(t) = £n=o M n 
with the boundary conditions /3(0) = f3(tf/2) = /3(tf) = 
-tt/2, /3(if) = -n/(2U), and /3(0) = n/(2tf), we can 
construct A and f2# [20[ . These two functions are shown 
in Fig. H for ^ = 0.2 ^s (Q R = 2X/h and A = -2Z/H). 
For the same process time tf we also plot in Fig. [5] JT (t) 
and Zo(t) for a Landau-Zener protocol in which the Rabi 
frequency is slightly larger than the maximun required 
for the invariant-based protocol: Qoj z = 30 MHz. As ex- 
plained in the previous appendix, an unreasonably high 
laser intensity would be required to make it adiabatic 
while satisfying the bare-state condition at the edges, 
and ilo,iz = 30 MHz is still too small to satisfy Eq. (|A3[) . 
This is evident in the failure to invert the population, see 
Fig. El We use a = -2800 MHz 2 to have the bare states 
as eigenvectors at the time edges which implies a rather- 
large detuning. Fig. [5] also depicts the Y(t) component 
of Hq 1 ^ and the X(t) and Z(t) components of for 
tf = 0.2 /zs. With these parameters these Hamiltonains 
provide shortcuts to adiabaticity, see Fig. |H1 but they use 
very high dctunings compared to those of the invariant- 
based protocol. This example does not mean, however, 
that invariant-based engineering is systematically more 
efficient. Invariant-based engineering and the counter- 
diabatic approach provide families of protocols that de- 
pend on the chosen interpolating auxiliary functions in 
the first case and on the reference Hamiltonian Hq in the 
second. Their potential equivalence was studied in [2fjj . 
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